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No snow-plough mechanism during the rapid hardening of 
supermassive black hole binaries 
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ABSTRACT 

We present two-dimensional hydrodynamical simulations of the tidal interaction between 
a supermassive black hole binary with moderate mass ratio, and the fossil gas disc where it 
is embedded. Our study extends previous one-dimensional height-integrated disc models, 
which predicted that the density of the gas disc between the primary and the secondary 
black holes should rise significantly during the ultimate stages of the binary's hardening 
driven by the gravitational radiation torque. This snow-plough mechanism, as we call it, 
would lead to an increase in the bolometric luminosity of the system prior to the binary 
merger, which could be detected in conjunction with the gravitational wave signal. We 
argue here that the snow-plough mechanism is unlikely to occur. In two-dimensions, when 
the binary's hardening timescale driven by gravitational radiation becomes shorter than 
the disc's viscous drift timescale, fiuid elements in the inner disc get funneled to the outer 
disc through horseshoe trajectories with respect to the secondary. Mass leakage across the 
secondary's gap is thus found to be effective and, as a result, the predicted accretion disc 
luminosity will remain at roughly the same level prior to merger. 
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1 INTRODUCTION 

The merger of two galaxies results in the formation of a newly 
assembled gala xy containing a superm assive black hole binary 
in its nucleus (|Begelman et al.lll98d ). Models of hierarchical 
galaxy formation predict that the most common supermas- 
sive black hole binaries arise from minor galactic mergers, 
and have unequal-mass black holes (jLacev fc Cole 19931 ). The 
[tidal interaction between the binary black hole, and the disc 
of gas and stars it is embedded in, shrinks the binary's or- 
bit. The hardening efli ciency has a complex d ependence on 
the binary's mass ratio (JKaza ntzidi s et al.''2003), the gas disc- 
to-binary mass r atio (jLodato et al.ii200gi"). and stellar dynam- 
ical interactions (|Milosavlievic fc Merrittil200ll ). Provided the 
binary hardens to sufficiently small separations, the angular 
momentum extracted from the binary's orbit becomes pro- 
gressively more dominated by gravitational radiation. When 
the gravitational wave torque prevails over the disc torque, the 
disc's accretion rate becomes much smaller than the binary's 
hardening rate. 

One-dimension al di sc models JArmitage fc NataraianI 

l2002l : iLodato et al]|2009l : IChang et al.fl2010l ) showed that the 
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density of the disc region located between the primary and the 
secondary black holes should significantly increase prior to co- 
alescence, as the binary's hardening is dominated by gravita- 
tional radiation. In this scenario, the secondary rapidly drains 
the inner disc onto the primary, acting like a snow-plough. 
The increase in the inner disc density would lead to a sud- 
den increase in the disc's bolometric luminosity in the few 
days prior to merger. This would constitute an electromag- 
netic precursor to a supermassive black hole merger, which 
could be identified in conjunction with the gravitational wave 
signal d etectable with the L aser Interferometer Space Antenna 
(LISA) l|Chang et al.ll2010l ). 

In this Letter, we show that this prediction does not hold 
in two-dimensional (2D) disc models. We demonstrate with 
the help of 2D hydrodynamical simulations that when the bi- 
nary's hardening is dominated by gravitational radiation, the 
inner disc is progressively funneled to the outer disc through 
horseshoe trajectories with respect to the secondary, with the 
consequence that there should be no significant increase in the 
disc luminosity prior to the merger. The physical problem we 
consider is described in § (2] and the numerical setup of our 
simulations is detailed in § [31 Our results of calculations are 
presented in § U followed by some concluding remarks in § [S] 
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2 PHYSICAL PROBLEM 

We examine the tidal interaction between a supermassive 
black hole binary with mass ratio q <^1 and the gaseous disc 
it is embedded in, neglecting for simplicity the presence of 
other stars. The torque responsible for the binary's hardening 
comprises the tidal torque Fdisc exerted by the gaseous disc, 
and the torque Few due to th e emission of gravitational wa ves 
(GW). The latter reads (e.g.. lArmitage fc NataraianlbOOsI ') 



32 G^ 

5 e'a^ 



G{Ah + M2) 
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Mi(Mi + M2)M| (1) 



for a circular binary, where Mi and Af2 ~ qM\ denote the 
mass of the primary and of the secondary, respectively, a is the 
binary's semi-major axis, c is the speed of light, and G is the 
gravitational constant. Since IFqwI strongly increases with de- 
creasing a, the binary's hardening is ultimately driven by the 
gravitational wave torque. The binary's hardening timescale 
TGw due to the GW torque is 
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Further denoting the Schwarzschild radius r^ 
Eq. ^ can be recast as 
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Toih qy/1 + q V ''g . 

where Torb is the orbital period at the binary's semi-major axis 
a. At the innermost stable circular orbit (ISCO), located at 
T'isco ~ 3rg for a non-rotating primary black hole, TGw/Torb ~ 
2.4, 25 and 250 for q = 10"S 10"^ and 10"^, respectively. 

In this study, the secondary to primary mass ratio is fixed 
to g = 2 X 10~^ , and this large value implies that the secondary 
is able to open a g ap around its orbit even in a relatively thic k 
and viscous disc (jLin Sz Papaloizoul 119861 : ICrida et al.l 120061 ') . 
It also makes the two-dimensional approximation for the disc 
applicable when discarding gas accretion onto the secondary 
D'Angelo et al.l 120051 1. We point out, as did iGould fc Rb3 
2000), that apart from the presence of the GW torque, the 
physical problem we consider shares a number of analogies 
with the orbital evolution of a massive gap-opening planet in 
a protoplanetary disc. 

We examine the evolution of the fluid elements in the in- 
ner disc close to the location of the inner separatrix of the 
secondary's horseshoe region. These fluid elements are on ap- 
proximately circulating streamlines with respect to the sec- 
ondary, with a relative (synodic) period Tsyn ~ half the horse- 
shoe libration period. The horseshoe libration period reads 
Tiib = Stto X (3ila;a)~^, where fl is the secondary's angular fre- 
quency, and Xs is the radial half- width of the horseshoe region . 
For gap-opening secondaries, Xb ~ 2Ru (JMasset et al.l 120061 ). 
where Ra = a(g/3)^''^ is the Hill radius of the secondary. 

During a synodic period, if the radial distance Sub by 
which the binary hardens becomes comparable to, or greater 
than the radial distance Saa by which fluid elements near the 
inner separatrix drift inward, these fluid elements embark on 
horseshoe trajectories and are funneled to the outer disc. The 
value of dua is essentially set by the viscous torque, and 6ab by 
Fdisc + Fgw ~ Fgw when the binary shrinkage is dominated 
by gravitational radiation. Thus, in order of magnitude, Sab ^ 
(5os when tgw becomes shorter than the viscous drift timescale 
~ 2a^ /Su, where u — ah^a^QK is the disc's turbulent viscosity 



{a denotes the alpha viscous parameter, h the disc's aspect 
ratio, and SIk the Keplerian angular velocity). Using Eq. Q, 
the typical binary separation below which Sab ^ Sas, often 
called decoupling radius, can thus be estimated as 
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where a and h are to be evaluated at the secondary's loca- 
tion. Detailed modeling of the disc's thermal balance yields a 
decoupli ng radius that is a few times la rger than that given in 
Eq. (gl) (|Milosavlievic fc Phinnevll2005l '). This illustrates that 
when the hardening of an unequal-mass binary black hole is 
dominated by gravitational radiation, a substantial funneling 
of the inner disc to the outer disc is possible beyond the ISCO 
location, which we confirm with 2D hydrodynamical simula- 
tions in § |4l 



3 PHYSICAL MODEL AND NUMERICAL 
SETUP 

We investigate the tidal interaction between a black hole 
binary and the gas disc it is embedded in, assuming the 
binary's hardening is dominated by the emission of grav- 
itational waves. For this purpose, 2D hydrodynamical 
simulations were carried out with the Fargo code (JMassed 
[2000 , http : //f argo . in2p3 . f r) . A cylindrical coordinate 
system {r, ip} centred onto the primary black hole is adopted. 
The disc extends from r = 1.2rg to r = 35.4rg (that is, 
r G [0.4 — 11.8]r-isco with assuming a non-rotating primary). 

Binary parameters — The mass of the primary and secondary 
black holes is Mi = 5 x 10^ Mp and M2 = qM i = W'^Mq, 
respectively, as in lArmitage fc NataraianI (|2002r ) . The initial 
semi-major axis of the binary is ao ~ ll.Srg « 5.7 x 10"** 
pc (or, ao ~ 4risco). The binary is assumed to be circular, 
and its orbital plane to be aligned with the disc. Each black 
hole is treated as a point mass potential. Relativistic effects 
are discarded for simplicity and to facilitate comparison with 
previous studies. Similarly, gas accretion onto the secondary 
is neglected. Also, to avoid a large accumulation of mass 
inside its Hill radius, the secondary's gravitational potential 
is smoothed over a softening length e — H{ao), where H is 
the pressure scale height (its value is specified below). 

Disc parameters — The disc is initially axisymmetric and 
rotates at the angular frequency Sl(r) about the primary. The 
expression for Q{r) assumes a radial equilibrium between 
the centrifugal acceleration, the radial acceleration due to 
the pressure gradient, and the gravitational acceleration 
due to the primary alone. The initial gas density is « 10^ 
gcm~'^ X (r/ao)~^'^, which corresponds to an initial mass 
Ri 6.5 X lO"** Ml. For simplicity, a locally isothermal equation 
of state is considered where the initial radial profile of the 
disc temperature, which we take proportional to r~^'^, 
remains stationary. The disc temperature is conveniently 
related to the disc's aspect ratio h — H/r, which we take to 
be h{r) = 0.08 x (r/ao)"'^'*. The disc turbulence is modeled 
with a constant kinematic viscosity v that corresponds to an 
alpha parameter a ~ 0.03 uniformly throughout the disc. 

Numerical setup — The grid has A'^,. — 280 cells along the 
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Figure 2. Time sequence of tiie gas surface density when the binary's hardening is driven by the GW torque. The x— axis shows the radial 
coordinate in units of the Schwarzschild radius rg (part of the grid is shown along this direction), and the y— axis shows the azimuth. The 
same color scale is adopted in all three panels. Time is in units of the binary's orbital period before the hardening stage. The binary's Hill 
radius is highlighted by a black circle, and a few streamlines in the frame rotating with the secondary arc overplotted by white curves. The 
location of the circular streamline of the inner disc closest to the secondary is shown by a red curve. 
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Figure 1. Time evolution of the binary's semi-major axis, a, driven 
by the GW torque alone (black curve) , and by the disc torque alone 
(red curve, both are results of simulations). Time is expressed in 
units of the binary's orbital period at its initial separation, and a 
is shown in units of r^, . 



radial direction, and it covers the full 2ti range in azimuth 
with Ntp = 780 cells. A logarithmic spacing is used along the 
radial direction to optimize the grid's resolution near the disc's 
inner edge. The secondary's Hill radius is resolved by about 
30 grid cells along each direction at all times. Standard zero- 
gradient outflow boundary conditions are applied at the grid's 
inner and outer edges. 



4 RESULTS 

Our sinmlations were performed in two steps. The secondary is 
first held on a fixed circular orbit about the binary's center of 
mass. Its mass is gradually increased over 100 orbits to avoid 
a violent disc relajcation following the introduction of the sec- 



ondary. The secondary progressively depletes a gap around its 
orbit, which attains a quasi steady-state density profile ~ 400 
orbits after the insertion of the secondary. Simulations were 
then restarted with the binary feeling the gravitational wave 
torque only. This second stage, which we refer to as the hard- 
ening stage in what follows, lasts < 410 orbits, as indicated by 
Eq. ^ . The disc torque on the secondary is discarded for sim- 
plicity, as IFdiscI ^ ITgwI fi'om the binary's initial separation. 
This point is clearly illustrated in Fig. [l] where we compare 
the time evolution of the binary's semi-major axis when only 
the GW torque is included (black curve), and when only the 
disc torque is included (red curve). 

Fig. [2] shows a time sequence of the gas surface density 
during the hardening stage. Time is expressed in orbital pe- 
riods at the binary's initial separation, ao- Streamlines in the 
frame rotating with the secondary are shown as white curves. 
The circulating streamline of the inner disc closest to the sec- 
ondary is depicted by a red curve. This time sequence clearly 
shows no density increase in the inner disc, but underlines in- 
stead the progressive replenishment of the gap. The gas den- 
sity inside the gap becomes more and more asymmetric as 
the binary's hardening gets faster, with more gas behind the 
secondary {ip < 0) than ahead of it. This asymmetry results 
from fluid elements of the inner disc embarking onto horseshoe 
streamlines and being funneled to the outer disc or to the gap. 
This gap asymmetry is reminiscent of the one occurring during 
type III runaway migration (a fast, gene rally inward migration 
regim e driven by the disc torque; see iMasset fc PapaloizoiJ 
l20(33l ). 

Fig. [3] displays the azimuthally-averaged disc density at 
different times prior to the binary's merger. It further illus- 
trates that the rapid shrinkage of the binary does not squeeze 



the inner disc in our 2D simulations, in contrast to 



previ- 



ous height-integrated ID disc models (lArmitage fc NataraianI 
l2002l : iLodato et al]|2009l : IChang et al.ll2010l ). 

To gain further insight into the time evolution of the mass 
in the inner disc, we have performed simulations in which a 
small patch of the disc inside the inner separatrix of the horse- 
shoe region is polluted initially with a passive contaminant 
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Figure 4. Time sequence showing the evolution of the specific concentration of the passive contaminant, C, defined in the text. The x— axis 
shows the radial coordinate in units of rg, and the j/— axis the azimuth. Contour levels are defined for each panel to underline the region 
occupied by the contaminant at each time. The binary's separation is indicated by a vertical dashed line. 
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Figure 3. Azimuthally-averaged surface density at different times 
expressed in units of the binary's initial orbital period. Radius is in 
units of the Schwarzschild radius rg. Density peaks correspond to 
gas accumulation inside the circum-secondary disc. 



that is advected with the flow. We denote by ij) the concentra- 
tion per unit area of the contaminant, and evolve the equation 



f + V-,*v, 



0. 



(5) 



Writing ip{r,(p) = C(r, </p) S, where C{r,ip) is referred to as 
the specific concentration of the contaminant, and using the 
continuity equation, Eq. ((5| can be recast as dtC ^w -V C = 0. 
The patch of contaminant is introduced 340 orbits after the 
beginning of the hardening stage, and the time evolution of its 
specific concentration, C, is displayed as a time sequence in 
Fig. |4l Time increases by « 7 orbits moving from left to right 
through the panels. The location of the secondary is shown 
by a vertical dashed line in each panel. In the left-hand panel, 
most of the contaminant is concentrated inside the inner sep- 
aratrix of the horseshoe region, located at r ~ 4.5rg. At this 
time, a (relatively) tiny fraction of the contaminant is being 
funneled to the outer disc, and ends up sliding along the outer 



separatrix of the horseshoe region at r « 8rg . The second panel 
highlights that more and more contaminant is being redirected 
to the outer disc. As the binary continues hardening, most of 
the contaminant now moving in the outer disc remains in the 
outer disc, only a small fraction of it undergoes inward horse- 
shoe U-turns back to the inner disc. In the third panel, most 
of the contaminant is now concentrated in the outer disc. The 
radial profile of the contaminant gets progressively thicker in 
the outer disc, because the outer separatrix of the horseshoe 
region is moving inwards faster than the contaminant may 
drift inward viscously. This is even more clear in the fourth 
panel of Fig. [l] where the contaminant in the outer disc is left 
well behind the secondary. We also mention that in addition 
to the material passing to the outer disc, we never observe the 
contaminant being pushed inward: the region below r ~ 3.5rg 
remains dark. 

The time evolution of the gas surface density in Figs. |2] 
and |3l and of the contaminant's specific concentration in 
Fig. 131 suggest that most of the mass in the inner disc rel- 
ative to the secondary is transferred to the outer disc through 
horseshoe streamlines. To help quantify this funneling mecha- 
nism, we display in Fig. [5]the time evolution of the mass in the 
inner disc, in units of the binary mass, during the last ~ 70 
orbits before the binary's merger. The increase in the mass 
of the outer disc during the same time span (also in units of 
the binary mass) is overplotted by a red curve. The difference 
between both curves accounts for the mass lost from the com- 
putational grid, which is shown by a dashed curve in Fig. [S] 
These curves clearly show that the fraction of the inner disc's 
mass funneled to the outer disc largely exceeds that drained 
onto the primary. 



5 CONCLUDING REMARKS 

The minor merger of two galaxies leads to the formation of a 
supermassive black hole binary with unequal mass ratio em- 
bedded in a gaseous disc. The tidal interaction between the 
binary and the disc hardens the binary's orbit. When the bi- 
nary gets sufficiently tight, emission of gravitational waves be- 
comes the main source of angular momentum extraction from 
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Figure 5. Mass of the inner disc during the ultimate stages of 
the binary's hardening, before the secondary leaves the computa- 
tional grid (black curve). The increase in the mass of the outer 
disc during the same time interval, namely AMouter = A-/outcr(t) — 
A^outerCSSTTo^jj), with Mouter the outer disc's mass, is displayed by 
a red curve. The mass lost from the computational grid during the 
same time interval is shown by a dashed curve. Masses are in units 
of the binary mass, and all curves have been smoothed over 5 orbits. 



iMacFadven fc Milosavlievid l2008l : ICuadra et ai] 120091 ') , mak- 
ing our funneling mechanism not applicable. Also, different 
values of the disc aspect ratio and viscosity would change the 
structure of the gap opened by the secondary (width, depth). 
Partial gap-opening could occur in very thick discs, and a 
three-dimensional disc modeling would be valuable. In some 
cases, the disc-secondary interaction can make the outer 
disc eccentric, thereby driving the bin ary's eccentricity (e.g., 
IPapaloizou et all 1200 ll : lOgilvie fc Lubo w 2003), although the 
latter should be partially damped d uring the fast inspiral 
drive n by gravitational radiation (|Armitage fc NataraianI 
|2005|). Still it remains to be clarified how the funneling 
mechanism would operate in the presence of an eccentric 
binary. 
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the binary's orbit, and causes further rapid shrinkage until 
coalescence takes place. 

We have focused in this Letter on the evolution of the 
disc region located between the primary and the secondary 
black holes (the inner disc), when the binary's hardening is 
dominated by the emission of gravitational waves. With the 
help of 2D hydrodynamical simulations, we have shown that 
the rapid hardening of the binary does not lead to a signifi- 
cant squeezing of the inner disc. The latter is redirected in- 
stead toward the disc region beyond the secondary's orbit (the 
outer disc) through horseshoe streamlines. When the binary's 
hardening timescale driven by gravitational radiation becomes 
shorter than the disc's viscous drift timescale, fiuid elements 
in the inner disc embark on horseshoe trajectories with respect 
to the secondary, and are progressively funneled to the outer 
disc. The funneling of the inner disc toward the outer disc 
implies that, in contrast to the predictions of ID disc mod- 
els, the accretion rate onto the primary is not dramatically 
increased just prior to merger, and, as a result, the disc emis- 
sion before the binary merger should remain at about the same 
level. After the merger, an electromagnetic afterglow could be 
detected as the disc gets accreted by the merged black hole 
(Milosavljevic & Phinncy 2005). 

The physical model we have considered is a straight- 
forward 2D extensio n of t he model c o nsider ed in 
Armitage fc NataraianI (|2002l ). iLodato et all l|2009l ). and 
Chang et al.l ( 2010f ). and as was already pointed out by 
these authors, this model has a number of simplifying 
assumptions. We have considered an intermediate binary 
mass ratio (g = 2 x 10"'^). A smaller mass ratio would 
decrease the binary's separation below which a significant 
funneling occurs (which we have checked with additional 
simulations, not reported here). A mass ratio closer to unity 
would initially lead to a rapid depletion of the inner disc, and 
the binary would be surrounded by a circumbinary disc (e.g., 
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